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Abstract. We present the first attempt of the systematical investigation about the effects of shell 
correction energy for dynamical processes, which include fusion, fusion-fission and quasi-fission 
processes. In the superheavy mass region, for the fusion process, the shell correction energy plays a 
very important role and enhances the fusion probability, when the colliding partner has strong shell 
structure. By analyzing the trajectory in the three-dimensional coordinate space with a Langevin 
equation, we reveal the mechanism of the enhancement of the fusion probability caused by shell 
effects. 

Keywords: Superheavy elements, fluctuation-dissipation dynamics, fusion-fission process, quasi- 
fission process 
PACS: 24.60.Ky, 27.90.+b 

INTRODUCTION 

In the heavy ion fusion reaction, with increasing atomic numbers of the target and 
the projectile, it becomes more difficult to make a compound nucleus due to the strong 
Coulomb repulsion force and strong dissipation force. In the superheavy mass region, 
this difficulty is more remarkable. Though the mechanism of fusion-fission reaction 
in the heavy mass-region is not clear, generally we recognize the existence of fusion 
hindrance. 

The fusion hindrance is mainly caused by the macroscopic properties of the colliding 
partner. On the contrary, the fusion is enhanced by the shell structure of the nuclei 
HI 0]. In the experiments, especially cold fusion reaction, these advantages are used 
to synthesize the superheavy nuclei 0,|4|]. To understand the fusion mechanism clearly, 
it is better to treat separately the fusion hindrance and the fusion enhancement, that is to 
say, the macroscopic aspect and the microscopic one. 

In our previous study ||5|], we have discussed the fusion hindrance in the superheavy 
mass region. We presented the origin of the fusion hindrance systematically by the 
trajectory calculation. In the present paper, as fusion enhancement, we focus on the 
influence of shell effects. It is known that the nuclear structure of the projectile-target 
combinations correlates to the touching probability |0, 0], but it influences also the 
dynamics from the touching point to the compound nucleus. We investigate precisely 
how the trajectory behavior is influenced by the shell correction energy. 
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FIGURE 1. Potential energy of the liquid drop model Vld (a) and with shell correction energy Vld 4- 
(b) for 292 1 14 in the z — (X space (5 = 0). The calculation is done by a two-center shell model code 



MODEL 

Using the same procedure as described in reference JH], we investigate the dynamical 
processes. The fluctuation-dissipation model with the Langevin equation is employed. 
We adopt the three-dimensional nuclear deformation space given by two-center parame- 
terization IE [HI. The three collective parameters involved in the Langevin equation are 
as follows: zo (distance between two potential centers), 8 (deformation of fragments) 
and a (mass asymmetry of the colliding nuclei); a = (A\ —A<i)/(A\ +Aj), where A\ 
and A2 denote the mass numbers of the target and the projectile, respectively. 

The multidimensional Langevin equation is given as 

dqi / _j\ 

— — = [m )-Pi, 
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where a summation over repeated indices is assumed, qi denotes the deformation co- 
ordinate, pi is the conjugate momentum of qt. V is the potential energy, and m, 7 and 
Yij are the shape-dependent collective inertia parameter and dissipation tensor, respec- 
tively. A hydrodynamical inertia tensor is adopted in the Werner- Wheeler approximation 
for the velocity field, and the wall-and-window one-body dissipation is adopted for the 
dissipation tensor. Details are explained in reference 




FIGURE 2. Sample of the trajectory projected onto the z — a plane at 8 = 0.0 (a) and the z — 8 plane at 
a = 0.0 (b) of Vld + Esheii f° r 284 1 14. The light yellow and white lines denote the trajectories which start 
at a = 0.46 and 0.24 at E* = 20 MeV, respectively. Symbols are given in the text, (c) Fusion probability 
with initial value a = 0.46 and a = 0.24, which are denoted by the red and blue lines, respectively. 



RESULTS 

Effect of the Cold Fusion Valleys in Fusion Process 

Figure [T] shows the potential energy surface of the liquid drop model (a) and with 
shell correction energy (b) for 292 1 14 in the z — OC space (5 = 0), which is calculated by 
the two-center shell model code ifTILUl. When we consider the shell correction energy, 



we can see the pronounced valleys which lead to the compound nucleus. The valleys 
are called 'cold fusion valleys' B 1 3L 1 1 411 . It is said that these valleys enhance the fusion 
probability. We discuss the effect of the cold fusion valleys in the dynamical process 
using the trajectory calculation. 

As our first attempt, in the dynamical calculation, we employ the potential energy 
of the liquid drop model with the full shell correction energy, which corresponds to 
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FIGURE 3. Potential energy surface V L d + E® hell in the z - a space at (a) 8 = 0.22 and (b) 0.40 for 
284 1 14, with a sample trajectory. 



the potential energy surface at the nuclear temperature T = MeV. It is represented 
by Vid + Vsh(1iT = 0), or V^o We compare the fusion probability with the 

different injection points, where the shell correction energy has remarkably negative 
and positive values like a = 0.46 and 0.24, that correspond to the reaction 76 Ge + 208 Pb 
and 108 Ru + 176 Yb, respectively. 

Figure 2(a) shows the sample trajectories which are projected onto the z — OC (8 = 0) 
plane of 284 1 14. The contact point is marked by (+). The light yellow line denotes the 
trajectory with the starting point a = 0.46 at the incident energy which corresponds to 
the excitation energy of the compound nucleus E*=20 MeV. This trajectory looks to 
move along the valley till z ~ 1.0. Then it enters the region near z ~ 0.0. On the other 
hand, the trajectory with the starting point a = 0.24 goes till z ~ 1.0 without changing 
the mass asymmetry parameter a, which is denoted by the white line. It seems that 
the trajectory overcomes the mountain located at z = 1.4, a = 0.24. We project these 
trajectories onto the z — 8 plane at a = in Fig. 2(b). The trajectory with the starting 
point a = 0.46 can enter around the compact shape region, but it is not trapped by the 
pocket around the ground state. The trajectory with the starting point a = 0.24 moves 
quickly in the +8 direction and goes to the fission region. We can say, the former is a 
deep quasi-fission process (DQF) and the latter is a quasi-fission process (QF) HE1]. 

Figure 2(c) shows the fusion probability with the initial value a = 0.46 and a = 0.24, 
which are denoted by the red and blue lines, respectively. The potential energy at the 
contact point with a = 0.46 is lower than that with a = 0.24, due to the shell correction 
energy. Therefore, at low excitation energy, the former fusion probability is larger than 
the latter one, because in the former case the available kinetic energy at the contact point 
is larger than in the latter case. The arrow denotes the Coulomb barrier ifl^l . 




FIGURE 4. Sample of the trajectory projected onto the z — 8 plane at a which corresponds to the 
turning point, (a) Vld case an d (b) VLD+E® hell case for 256 No. Vu)+E® he[l with a = 0.0 is denoted in (d). 
The evolution of nuclear shapes is presented in (b). 



Evolution of the Cold Fusion Valleys in the Dynamical Process 

Usually the cold fusion valleys are discussed using the potential energy surface on 
the z — OC plane with 5 = 0, like as Fig. 1(b). However, in the dynamical process, the 
trajectory moves in the large +8 direction. We should discuss the cold fusion valleys 
with the dynamical evolution of the 8 parameter. For the 5 = case, which corresponds 
to Fig. 1(b), we can see the remarkable the cold fusion valleys. With changing 8 value 
in the dynamical process, the cold fusion valleys also changes, which is shown in Fig. 3. 
When the 8 changes about 0.4, the cold fusion valleys disappear. 



The Role of Shell Effects in the Fusion Process 



With the shell correction energy, we discuss another important effect to enhance 
the fusion probability. In the discussion on fusion hindrance [5], we indicated that the 
turning point is important. Fig. 4(a) shows the potential energy surface of Vld for Z=102 
at the turning point on the z — 8 plane, that is to say, a corresponds to the value at the 
turning point. The mean trajectory is denoted by the white line. The gray line denotes the 
trajectory with taking into account the fluctuation. At the turning point (z ~ 0.3, 8 ~ 0.5), 
the trajectory moves to the fission direction due to the potential landscape. 

On the other hand, when we take into account the shell correction energy in Fig. 4(c), 
the temporary pocket appears at the turning point (indicated by A). At a = in Fig. 4(d), 
the pocket (indicated by B) corresponds to the ground state. When we take into account 
the shell correction energy, the trajectory is trapped in the pocket A at the turning point, 
and is blocked from going to the fission area II 1 611 . During the stay in the pocket A, 
the mass asymmetry is relaxed. With appearing of the large pocket B at a ~ 0, the 
trajectory moves to the pocket B. The pocket A helps the trajectory to enter the region 
corresponding to the compound nucleus. Fig. 4(b) shows the time evolution of the 
nuclear shape. The horizontal axis denotes the 8. The trajectory moves in the large 8 
direction quickly, then it is trapped in the pocket A. During the stay in the pocket A, the 
mass asymmetry changes. With approaching a = the trajectory moves into the pocket 
B. It corresponds to the compound nucleus. 

The temporary pocket that appeared due to the shell correction energy enhances the 
fusion probability. As the further study, we take into account the temperature dependence 
of shell correction energy. 

The author is grateful to Professor M. Ohta for his helpful suggestions and valuable 
discussion throughout the present work. 
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